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Abstract 

Using the Kondo lattice model with classical spins in infinite dimension, mag- 
netic phase transition in the perovskite-type 3d transition-metal oxide Lai_ x Sr x Mn03 
is theoretically studied. On the Bethe lattice, the self-consistency equations are 
solved exactly. Curie temperatures at the region of double-exchange ferromagnetism 
0.1 ^ x & 0.25 as well as the Neel temperature at x = are well reproduced quan- 
titatively. Pressure effect on the Curie temperature is also discussed. 
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After the discovery of high T c oxides, the physics of strongly correlated 3c? transition- 
metal oxide compounds have been revisited extensively. One of such materials is the 
perovskite-type manganese oxides Lai_ x Sr x Mn03. The most prominent feature of the 
system is the giant magnetoresistance (GMR) with negative sign. The system is a double- 
exchange ferromagnelfi!) at x ^ 0.1, while the antiferromagnetic phase with spin canting 
is observed at x ~ 0. For the detail of the phase diagram and the GMR in Lai-xSr^MnO^ 
and related materials, the readers are referred to e.g. refs. [3] and [|, and those cited therein. 

In Lai_ x Sr x .Mn03, 3c? electrons are considered to form both localized spins with S = 
3/2 in the t2 g orbitals and itinerant electrons in the e g orbitals which are coupled each 
other by Hund's ferromagnetic interaction. The bandwidth of the itinerant electron is 
estimated to be W ~ leV from the recent band calculation,^ while the Hund's coupling 
is considered to be larger than the bandwidth. As a model Hamiltonian of this system, 
the Kondo lattice model (KLM) with S = 3/2 and the ferromagnetic exchange coupling 
has been proposed.!^ 

From the strong coupling limit of the above model, Curie temperature T c has been 
studied.i'i) Using a mean-field type treatment, Curie temperature is estimated to be 
T c ~ W, where W is the bandwidth of the conduction electron. This result is roughly 
understood as follows; the transition temperature is determined from the competition 
between the gain of kinetic energy in the ferromagnetic state and the gain of entropy in 
the paramagnetic state. 

Quantitatively, however, the mean-field approach fails to explain the experimental 
value of T c , e.g. T c ~ 300K at x ~ 0.2 which is a few orders of magnitude lower than 
the theoretical estimate. The discrepancy seems to be due to the fact that the mean-field 
treatment poorly describes the electronic state in the paramagnetic phase so that kinetic 
energy is not estimated properly. In the paramagnetic phase, thermal spin fluctuations 
have to be taken into account appropriately because the system is in a strongly correlated 
regime. Recently, Millis et.al. has made an analysis on T c based on the calculation of 
the spin stiffness at the ground stated Their result also overestimates T c . Since the 
changes in the electronic state from the ferromagnetic ground state to the paramagnetic 
state should be nontrivial and drastic, ground state properties may not reproduce the 
transition temperature. Because the GMR is most prominently observed in the vicinity 
of the Curie temperature, a theoretical approach which is able to predict T c with accuracy 
is demanded also from the standpoint of a material designing and applications. 
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The author has shown in his previous wor 



kPi'S that the KLM in infinite dimensional 



limit D = oo and infinite high-spin limit S = oo reproduces the transport properties of 
Lai_ x Sr x Mn03 quantitatively. Within these limits, Green's functions are obtained exactly 
even in the paramagnetic state. Therefore, a proper treatment is possibly performed to 
obtain magnetic instabilities in the paramagnetic phase. In this paper, we study the 
magnetic phase transition of the above model. One of the aims is to examine, within this 
model, whether it is possible to explain the magnitude of the Curie temperature as well 
as its doping dependence in a consistent way. 

Thus we investigate the KLM at D = oo and S = oo. The Hamiltonian is described 

as 



where = (m^, rrii v , mf) and \m\ 2 = 1. We consider the Bethe lattice in the infinite 
coordination number limit. From the nature of the Bethe lattice that lattice points are 
divided into two sublattices, it is possible to study magnetic phases with ferromagnetic 
and antiferromagnetic order parameters. The density of states (DOS) is given in the 



semi-circular form N (e) = (2/irW)Jl — (e/W) 2 , where W is the bandwidth. Here we 
see another advantage of the Bethe lattice that the DOS has similar properties to that 
of the D = 3 systems, namely the existence of a band edge and the shape of the DOS 
at the edge. The above properties may be important in the case of making comparisons 
between the theoretical results and the experimental data in the low hole concentration 
region. 

The infinite-dimensional model is mapped to that of the corresponding single-cell 
model. The partition function is written in a form 



where integration is performed with respect to the orientation of local spins. Using the 
Weiss fields Go a (iu n ) which describe the dynamical motions of electrons for a = A, B 
sublattice sites, the fermion trace is given by 




(1) 



<ij>,a 





(2) 



Z { (rhA,rh B ) = J\ 4exp J^logdet 



a=A,B \ n 



(G 0a + Jm a a) 



) 



(3) 
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The Green's function is calculated exactly as 



Z { (rn,A,rn B ) i^-x. 



G a (iu n ) = J J dn A dn B IV | D > (G^(ko n ) + Jm a a) . (4) 
On the Bethe lattice with two-sublattice symmetry, the mapping relation gives0) 

G 0A ~ l {\uj n ) = iu n + fi- G B (iuj n )/A, 

GQ B ~ l (iu n ) = iu n + fi - G A (ko n )/i. (5) 

Self-consistency equations (0)-(H) are mainly solved in a numerical way. At J/W ^> 1, 
Green's functions near the magnetic transition temperature are also studied analytically 
Sublattice magnetization of the local spins (rh a ) is obtained from 

/ - \ f f An An - Z f (m A ,m B ) 

K) = J J dil A dil B m a . (6) 

Then, magnetic transition temperatures are determined as a function of Hund's coupling J 
and chemical potential /x in the unit of W. The carrier electron number for Lai_ x Sr x Mn03 
is nominally considered to be n = 1 — x. Hereafter we use the hole picture so that the 
hole concentration is expressed by (x) — 1 — (n). In this paper, we restrict ourselves to 
< (x) ^ 0.5, i.e. from half-filling to quarter-filling. 

In Fig. [I], we show T c /W as a function of J/W at various filling. Here, for simplicity 
of the calculation, chemical potentials are systematically chosen to be //(J) = —J + 5jj, 
with Sfi/W = 0, 0.25 and 0.33 so that the carrier numbers at T c become (x) ~ 0.5, 0.3 
and 0.2, respectively, with errors ^ ±0.03 at J/W > 4. The curves in the figure are the 
results of the (J/W) -1 expansion at J/W 3> 1 in the form 

TJW = f [l - JiJ/W)- 1 ] . (7) 

We see T c cx W at J/W — > oo. Dimensionless constants are given by T ~ 0.044, 0.039 
and 0.035 while J ~ 0.49, 0.89 and 1.12 for (x) ~ 0.5, 0.3 and 0.2, respectively. The 
fitting curves seems to reproduce results also at J/W ~ 1 surprisingly well. The quantity 
T is the Curie temperature in the scale of W at J/W —>■ oo, while J may be interpreted 
as the magnitude of the instability of the ferromagnetic order upon decreasement of J/W. 
From T ~ 0.04 and W ~ 10 4 K, we see that the Curie temperature in Lai_ x Sr x Mn03 is 
roughly explained within this model in the strong coupling limit. The tendency that T 
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(J) increases (decreases) as holes are doped is easily understood because the stability of 
the ferromagnetic order should be enhanced by hole doping. 

Now, we precisely compare the doping dependence of T c with the experimental data. 
Calculations are performed at J/W — 4 which is the value that explains the magnetore- 
sistance phenomena successfully,!!) as well as at J/W = oo. In Fig. |, we show T c as a 
function of x together with the experimental data from ref. |j. Here, the bandwidth is 
scaled so that T c at x = 0.15 reproduces the experimental result; we then have W = 1.05eV 
for J/W = 4, which is a moderate value. Thus we see from the above result that the 
experimental data are quantitatively reproduced very well at x ^ 0.25. Discrepancy on 
the doping dependence of T c is observed at x ^ 0.3, which may be interpreted as a sign 
of a crossover of the system from the strong coupling regime to a weak coupling limit. 

Next, we calculate the pressure effect on T c . In Lai_ x Sr x Mn03 at the ferromagnetic 
metal region x ^ 0.15, increase of T c is observed under pressure.0-* The pressure coef- 
ficient dlogT c /dP is positive, and it decreases systematically as the hole concentration 
is increased. For the theoretical treatment, we make following assumptions; the band- 
width increases as pressure is applied due to the increase of overlaps between neighboring 
orbitals, while the intra-atomic Hund's coupling is not affected. Then, we have 

dlogT c = dlogiy<91ogT c 
dP dP dlogW' [ ' 

where d\ogW/dP > from the assumption. Using eq. (|7|) at J ^> JW, we have 

1 - J{J/W)-\ (9) 



^logv; . , Tl r/Tr .._, 



dhgW 

so that a positive pressure coefficient is derived in the strong coupling region. We have 
shown above that J decreases as x is increased. Then, the pressure coefficient of T c should 
increase by hole doping if J/W is fixed, which is contrary to the experimental result. In 
order to explain the experiment within the present model and the assumptions mentioned 
above, we must have the effective decrease in J/W as the hole concentration is increased. 
This is another implication that the system undergoes the crossover to the weak coupling 
region by hole doping. 

In Fig. |^, the Neel temperature Tn at x = is shown. At J ^> W, the expansion with 
respect to (J/W)^ 1 gives an analytical formula 

T N /W = (1/24) • (J/W)-\ (10) 
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The Neel temperature in LaMnOs, Tn — 140K,eP is also well explained, e.g. if we set 
J = 4W and W ~ 1.2eV. Upon hole doping, the model in the low temperature regime 
shows a first order transition from the Neel state directly to the ferromagnetic state. 
Canted Neel state which is observed in Lai_ x Sr x Mn03 at x ^ 0.1 is not found within this 
model at J/W ^ 2. It must be noted that the present approach does not take into account 
the magnetic phase with incommensurate wave numbers such as spiral states or conical 
states. In order to study the magnetic structure at x ^ 0.1, it seems to be necessary 
to treat the model in a realistic lattice structure. Effects from the orbital degeneracy in 
e g bands and its Jahn- Teller splitting may also be essential to explain the experimental 
results in that region. 

To summarize, we have calculated the KLM in the limit S = oo and D = oo on 
the Bethe lattice. Magnetic transition temperature is obtained exactly as a function of 
interaction strength and hole concentration. Comparison with the experimental data in 
Lax_2.Srj.MnO3 is made. At J/W ~ 4 and W ~ leV, the model consistently explains the 
Neel temperature at x = and Curie temperatures at 0.1 ^ x ^ 0.25 in a quantitative 
way. 

The author would like to thank Y. Tokura and T. Arima for fruitful discussions and 
comments. The numerical calculation is partially performed on the FACOM VPP500 at 
the Supercomputer Center, Inst, for Solid State Phys., Univ. of Tokyo. 
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Figure Captions 

Fig. 1: Curie temperature T c /W as a function of J/W. 

Fig. 2: Curie temperature T c at J/W = 4 and J/W = oo as a function of x. Exper- 
imental data in Lai-xSr^MnOs are from ref. |j. Temperature is scaled in the 
unit of Kelvin. 

Fig. 3: Neel temperature at x — 0. The line in the figure is the result from (J/W) -1 
expansion. 
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